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Abstract 

An interesting type of spin chain has appeared in the context of the 
planar AdS/CFT correspondence: It is based on an integrable nearest- 
neighbor spin chain, and it is perturbatively deformed by long-range 
interactions which apparently preserve the integrable structure. Sim- 
ilar models can be constructed by demanding the existence of merely 
one conserved local charge. Although the latter is not a sufficient inte- 
grability condition in general, the models often display convincing signs 
of full integrability. 

Here we consider a class of long-range spin chains with spins transform- 
ing in the fundamental representation of qI{N). For the most general 
such model with one conserved local charge we construct a conserved 
Yangian generator and show that it obeys the Serre relations. We thus 
provide a formal proof of integrability for this class of models. 



*Work performed at Joseph Henry Laboratories, Princeton University, Princeton, NJ 08544, USA 



1 Introduction 



Among integrable spin chains the models with pairwise interactions of spins at adjacent 
sites are understood best. Their integrable structure is based on an R- matrix from 
which many interesting quantities and their relations can be derived. For instance, the 
R-matrix gives rise to commuting transfer matrices from which a Hamiltonian and a set 
of commuting conserved local charges follows. The existence of the latter is one way to 
define the integrability property of a Hamiltonian. 

Let us consider the simplest type of model with manifest Lie algebra symmetry. The 
R-matrix of such models typically is rational and it is well known that it possesses an 
additional hidden symmetry: the Yangian. The Yangian is a Hopf algebra which extends 
the Lie symmetry to an infinite-dimensional algebra (see [2] for reviews and further 
references) . Formally, the Yangian is a deformation of the universal enveloping algebra 
for the (half) loop algebra of the Lie symmetry. The Lie generators, J ^ reside at level-0 
of the loop algebra and we shall call the level- 1 generators, the Yangian generators. 
Repeated commutators of the latter supply us with all the higher levels of the loop 
algebra and therefore there is no need to consider them explicitly. Merely, the Yangian 
algebra imposes a constraint on triple products of generators: the Serre relations. 

Just as the Lie generators, the Yangian generators have a representation on the spin 
chain. Whereas the Lie generators act locally along the chain, 

^^ = Y.,^t (1-1) 

the Yangian generators act as bi-local products of Lie generator insertions 

= E.<, l^iaJk^^Jr- (1-2) 

Here Jj^ is the Lie generator J'^ acting on site k of the chain and F^q are the structure 
constants of the Lie algebra. The Yangian typically is not an exact symmetry of the 
Hamiltonian, it merely commutes up to contributions at the ends of the spin chain. Al- 
ternatively one can say that the boundary conditions normally break Yangian symmetry. 
Nevertheless, for an infinite spin chain the boundary contributions are dropped and the 
Yangian becomes an exact symmetry also of the Hamiltonian. Yangian symmetry (up 
to boundary terms) is thus another way to define the integrability of a Hamiltonian. 

The integrable structures for spin chains with long-range interactions are less well 
understood. Among these, the best known class of models are the Haldane-Shastry 
chains [3]. The Hamiltonian of these models consists of interactions between two spins at 
a distance, and it is actually exactly invariant under a Yangian algebra. Less clear are the 
Inozemtsev spin chain models, which are a generalization of the Haldane-Shastry models, 
but whose Hamiltonian commutes with the Yangian only up to boundary contributions. 

A much wider class of models [IHTj was found in the context of integrability in the 
so-called AdS/CFT correspondence (string/gauge duality) [8], cf. [9ll3lll0] and [H] for 
reviews, and in the plane- wave matrix model [I2], cf. [l3l[Tl]. Here the Hamiltonian 
consists of interactions of arbitrarily many neighboring spins. An ordering principle for 
the interactions is provided by a global coupling constant A which is assumed to be small. 
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The Hamiltonian is given as a series expansion in this couphng 

oo 

H(A) = ^A'h('). (1.3) 

The leading-order Hamiltonian TC^^^ is of nearest-neighbor type and the range of inter- 
actions increases by one site per order in A: Ti^^-* consists of interactions among i + 2 
sites. These models were called long-range spin chains because for finite values of the 
coupling the range of interactions would be infinite. On the other hand, one might also 
call them short-range models because up to a given perturbative order (as much as can 
be constructed explicitly in practice) the interactions are local and the range is bounded. 
Note that the Inozemtsev models [15] are effectively within this class [16j, but for our 
models the interactions are more general and can be between more than two spins at the 
same time. Also the Hubbard model can effectively be described by a short-range spin 
chain in a certain limit [17J. 

The investigation of short-range integrable spin chains was initiated in [4j| and contin- 
ued in [5],[T8l[6]. An obvious complication for this class of models is that there exists no 
general framework to construct and to prove the integrable structures. Unlike the case 
of nearest-neighbor models one cannot rely on an R-matrix because it could generate 
only interactions between two spins at a time and not between arbitrarily many. Instead 
one had to explicitly construct some commuting local charges and rely on good faith 
that infinitely many further charges exist. Intriguingly it turned out that in all cases 
considered the existence of merely one commuting local charge was sufficient to ensure 
the existence of as many further charges as could be constructed in practice. 

The next immediate challenge was to determine the spectrum of short-range spin 
chains. Based on the special case of Inozemtsev models p^llG]. a proposal for the 
Bethe equations for a special short-range model was made in [19]. Furthermore, certain 
deformations of the Bethe equations [20] could be identified with certain deformations 
of the Hamiltonian [3] . These deformations are important for an accurate description of 
the AdS/CFT integrable model [TJEI] and a sector of the plane- wave matrix model |14j . 
Finally, a thorough derivation of Bethe equations for short-range models based on the 
asymptotic coordinate Bethe ansatz [18] was proposed in [6]. These Bethe equations 
have the curious feature that they can be written down for arbitrary finite values of 
the coupling A while the Hamiltonian was only known at some perturbative order. One 
should note though that the Bethe equations are merely asymptotically correct and 
not exactly valid as for nearest-neighbor models. This means that the Bethe equations 
predict the spectrum of a chain of length L correctly only at the perturbative order 
i = L — 2. The currently most advanced study for short-range models with 0[(A^) 
symmetry and spins in the fundamental representation [S] constructed the most general 
third-order Hamiltonian and made a proposal for the complete moduli space at arbitrary 
order and the corresponding asymptotic Bethe equations. 

Nevertheless, the class of short-range integrable spin chains is still in its infancy. 
For example, we still have no proof of integrability beyond next-to-leading order, we 
merely have indications due to the construction of one or more conserved charges. Yan- 
gian symmetry would be one way to produce a formal proof of integrability at a given 
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perturbative order. It is conceivable that short-range Hamiltonians display such a sym- 
metry (up to boundary terms). Indeed there are indications that this is true in some 
cases • However, many of the models considered in these works are equiva- 

lent to Inozemtsev models at the given accuracy, and therefore we cannot tell if the most 
general deformation investigated in [6] still has Yangian symmetry. In fact, the most 
general models have scattering matrices [T81I25] which are not functions of the difference 
of some rapidity variables. This may or may not be an obstacle for Yangian symmetry 
which commonly leads to scattering matrices of difference form. 

Consequently, we would like to see if general short-range models do display Yangian 
symmetry. This might also teach us more about the algebraic structures underlying 
short-range integrability. In this work we shall focus on the models with qI{N) symmetry 
and spins in the fundamental representation. 

We start in Sec. [2] with a review of short-range spin chains, more specifically we 
review the most general perturbative deformation with gl{N) symmetry studied in [H] 
along with the asymptotic Bethe equations for such a spin chain. In Sec. [3] we review 
Yangian symmetry of spin chains which is used later to prove integrability. The Yangian 
for the Inozemtsev chain is discussed to motivate the deformations of the Yangian in 
our spin chain. In Sec. H] we discuss the construction of the Yangian. With all the 
machinery ready we proceed to the calculation. We compare three different methods of 
constructing the system; these methods differ only by the order in which the constraints 
are applied. We then discuss the implications our findings. In App. |A]we present some 
of the simplifications used in carrying out the calculation and in App. [B] we explicitly 
perform the procedure to 0{\) for concreteness. 

2 Short-Range Spin Chains 

Here we introduce our perturbatively integrable short-range spin chain models with gl{N) 
symmetry introduced in [1]. We merely review the setup and some results, further details 
of the calculation can be found in the literature [HE]- 



We consider a spin chain with spins transforming in the fundamental representation of 
Qi{N). Thus a spin chain state of length L is an element of (C^)®^. A basis for such 
states is given by 



with ak = 1, . . . , N. 

We write the Hamiltonian as a power series in the coupling constant A as in fll.3p . At 
leading order the Hamiltonian TC^^^ is a nearest-neighbor interaction as for conventional 
spin chain models. For each perturbative order the range of the interactions may grow 
by one site. Thus the order A'^ correction Ti.^'^^ consists of interactions of range k + 2. 

This class of models is quite general. As a special case it includes the Inozemt- 
sev hyperbolic chain [15] when the parameter is expanded around the point where the 
Hamiltonian becomes a nearest-neighbor interaction pTO]. Furthermore the short-range 



2.1 Setup 




(2.1) 
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spin chain which arises in the su(2) sector of planar A/" = 4 supersymmetric Yang-Mills 
theory [ll[6] is of this kind. 

The higher local conserved charges of the integrable model are denoted by Qr- They 
have a similar perturbation series as the Hamiltonian 

oo 

Qr = Y.^'^r'\ (2.2) 

fc=0 

but the range is increased by r — 2 units. In other words, Qi'^^ consists of interactions of 
range k + r. For an integrable system all the charges have to commute with each other 

[Qr, Qs] = 0. (2.3) 

The Hamiltonian is defined to be the lowest charge Ti := Q2- In this paper we will 
only consider the next higher conserved charge and denote it by Q := Q3 for simplicity. 
Consequently, the defining relation for the integrable system is [7i, Q] = 0. 

The Qi{N) generators will be denoted by J^^. They act on a spin chain as follows 

L 

J'fTIci, . . . ,cl) = J^^fc^ci, . . . ,Cfc_i,a,Cfc+i, . . . ,cl). (2.4) 

k=l 

The Hamiltonian and the charges are assumed to be invariant under Qi{N) 

[J,\n] = [J,\Qr]=Q. (2.5) 

2.2 Notation 

We use the symbol [ai, 02, . . . , ctn] to describe a local 0[(A^)-invariant interaction of range 
n. The numbers ai, 02, . . . , are distinct integers between 1 and n which describe the 
permutation vr with TT{k) = a^. The interaction sums over all sets of n adjacent sites of 
the spin chain and permutes them according to vr: 

). (2.6) 

k 

For example [1, 2, 3, . . . , n] (for any n) acts identically on the state and essentially counts 
the number of n adjacent sites within the chain: This number equals the length L on 
a periodic chain. The symbol [2, 1] sums over nearest-neighboring spins and permutes 
them. This is the essential contribution of the Hamiltonian of a q\-{N) nearest-neighbor 
spin chain. A general local operator can be visualized as in Fig. [H 

Note that there is an implicit sum in this notation since the position of the first leg 
is summed over all sites of the spin chain. With this implicit sum in mind we see that 
many of the operators are equivalent: 

[ai, a2, . . . , a„, n + 1] = [ai, 02, ... , a„] = [1, ai + 1, 02 + 1, . . . , a„ -M], (2.7) 

since the last or first leg, respectively, acts trivially. 
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Figure 1: Graphical representation of a local invariant interaction of range n. 
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Table 1: Free parameters. 



The interactions [ai, . . . , a„] are all invariant under 0t(iV) and they form a complete 
basis of invariant interactions of range n. However, due to the identification in (12. 7p a 
minimal basis for range n consists of the interaction of range m < n with all leading 
leg Oi 7^ 1 and trailing leg a„ 7^ n. It can be shown that there are n! — (n — 1)! + 1 
inequivalent such interactions of range n. 

2.3 Construction 

The most general integrable Hamiltonians with the above properties were constructed 
in [6] , let us review the results. At the leading order, the Hamiltonian and the charges will 
take the standard form: a nearest-neighbor integrable spin chain with qI{N) symmetry. 
At a fixed order in perturbation theory we will make the ansatz that the Hamiltonian 
Ti,^^'^ and the conserved charges Q^^^ are linear combinations of all the local invariant 
interactions of the aforementioned range. The linear combinations are given by a set of 
a priori undetermined coefficients. We then impose that Ti and Q commute at every 
order in perturbation theory 

k 

^[^W^g(fc-^)] = 0. (2.8) 

At any given order this system of equations is linear in the coefficients of Ti^*^-* and Q'-'^-*. 
It turns out that, conveniently, it suffices to solve for these coefficients only. The contri- 
butions at intermediate orders, < i < k, serve as inhomogeneities for the system. The 
solution to the system has several free coefficients a, 7, e; their counting is summarized 
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in Tab. [H These correspond to solutions S'H^''\ SQ^^^ of the homogeneous system 

[6n^^\ g(o)] + [^(0)^ ^g(fe)] ^ 0. (2.9) 

The main task addressed in [6J was to understand the origin of the various coeffi- 
cients. While most coefficients appear in both STi.^''^ and SQ^''\ some coefficients can 
be attributed to only one of them. These stem from solutions of the stronger equations 
[5'H^''\ Q^^^] = or ['H^'^\ SQ^'^^ = 0. Obviously the solutions are the leading-order con- 
served charges. In the full perturbative setup they correspond to the mixing of conserved 
charges 

oo oo 

n = ^2,cC + Y,l2,s Qs, Q = l3,C^ + Y.^3,s Qs, (2.10) 

s=2 s=3 

where Qs are some bare commuting charges and where C is the length operator C = [1]. 
Clearly the Hamiltonian H and the charge Q obey the same properties as the Qs provided 
that the coefficients start at the proper order 

^r,C = 0(A°), 7r,s = 0(A--^(°'^-^)). (2.11) 

After subtracting these degrees of freedom, the remaining coefficients are strictly 
shared between Ti, and Q. We will call them intrinsic coefficients of the integrable system 
and discuss them in the next subsection. 

The final results for the bare charges Q2, Q3 were already obtained and presented in 
the earlier work [6], however in a different notation. We present them in our notation 
up to second order in A at the end of the paper in Tab. [6|8l For the convenience of 
the unexperienced reader we perform this process explicitly in App. [B]up to C(A). The 
higher order calculation is done precisely the same way but with significantly more terms. 



2.4 Bet he Equations 

The intrinsic parameters a^, Pr,s, ^k,i in the operators Tab. [6)18] can be understood best by 
considering the Bethe equations. The asymptotic Bethe equations (taking the assump- 
tion of integrability for granted) were derived in [B]. They read 

~f r^l 11 : exp(2«0 Mi,fc,Mij ) M — . 

^1-1 i Ki Ki^i j 

^ _ -TT ue,k — Ue-ij — 2 TT- U£^k — ue,j + 1 -rr 'Ui,k — Ue+ij — 
jJl ULk - + I Uii^k - Uej - i |jL u^^^ _ . _^ 1 

Kn~2 i Kn-l 

^ _ TT '^n-l,k — ^n-2,j — 2 TT Un~l,k — ""n-lj + ^ ^2 ]^2) 

Un-l,k - Un-2,j + \ Mn-l,fe " Mn-l,i - % 
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We now explain the rapidity map x{u) and the dressing phase 6{u,v). 
The rapidity map x{u) is defined imphcitly through its inverse 

The inverse map from the w-plane to the x-plane has the following form 



, , u u 
An) = - + - 



u 



The parameters adX) are fixed uniquely by the components of afc(A). The coefficients 
OLiiX) govern the propagation of spin fiips in the ferromagnetic vacuum (magnons). 
Next we present the dressing phase 



oo oo 



9{u,v) = '^ ^ Pr,sW{qr{u) Qsiv) - qs{u) qr{v)) . (2.15) 



r=2 s=r+l 

These functions govern the scattering of two magnons. For fl2.15p we also define the 
elementary magnon charges as 

The solutions to the above asymptotic Bethe equations define the set of periodic 
eigenstates of the system. Finally, we have to specify the eigenvalues of the Hamiltonian 
and higher charges in terms of the rapidities fc. The eigenvalue of the shift operator 
exp(2P) is given by 

In particular, cyclic states obey the zero- momentum condition exp(zP) = 1. The eigen- 
values of the bare spin chain charges are determined by the formula 



Qr = J2lrM. (2.18) 



k=l 



In the perturbative scheme the functions a^(A) and /3r,s(A) are series in A starting at 
the orders 

a,(A) = 0(A^+i), Pr,sW = O(A^-i). (2.19) 

For the coefficients Q;fc(A) it implies that ao(A) = C(A) and ae>i{X) = C(A[^/^1+^). 

The coefficients e^^i do not appear in the Bethe equations and therefore they can 
have no impact on the spectrum. Instead they correspond to perturbative similarity 
transformations of all the operators 

oo 

Qr = rQrr-\ r = 1 + ^ A'^r^'^^ (2.20) 

k=l 
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Table 2: Intrinsic parameters of the short-range chain. 



where T'-*^-' is an arbitrary interaction of range k + 1 parametrized by ek,i- Note that 
contributions to T which are linear combinations of the commuting charges do not alter 
the charges. Thus for counting purposes one has to remove these trivial similarity trans- 
formations. The total counting of intrinsic coefficients is collected in Tab. [21 The table 
shows that up to 0{X^) all the intrinsic parameters can be accounted for and we have 
thus fully understood the moduli space of short-range integrable spin chain models with 
Qi{N) symmetry. 

In fact we merely made the integrability of the model plausible, but have not explicitly 
proved it. The two known signs of integrability for this model are: one conserved charge 
exists, and the construction of further charges (as far as it has been tested) has in no case 
led to additional constraints on the Hamiltonian. The spectrum derived from the Bethe 
equations (as far as it has been tested) agrees precisely with the diagonalization of the 
Hamiltonian, see [H]. In the next section we aim for a proof of perturbative integrability 
by means of Yangian symmetry. 

3 Yangian Symmetry 

For integrability we need one conserved quantity per degree of freedom which seems hard 
to show explicitly for the model introduced in the previous section. Instead we can show 
that the spin chain Hamiltonian commutes with a Yangian algebra (up to boundary 
terms). The algebra implies the existence of an infinite tower of conserved generators. 

3.1 Yangian Algebra 

We will deal with only the lowest two levels of this algebra: the lowest level contains the 
Lie symmetry generators JT, in our case those of g = gl(A^). We shall call generators at 
the next level the Yangian generators and we shall denote them by 3^. All the generators 
at higher levels can be constructed iteratively from commutators of y provided that 
the generators obey the Serre relations of the Yangian. Since both J', y commute with 7i, 
the higher generators will also commute and we have an infinite tower of conserved (non- 
local) charges. See for a derivation of the Serre relations and the Yangian algebra. In 
order to prove integrability of the spin chain Hamiltonian it thus suffices to show: 

• Conservation of a Lie-type symmetry q, i.e. [H, J'] = 0. 

• Conservation of a g-adjoint generator 3^, i.e. [H, y] = 0. 
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• The two sets of generators J ^ y must obey the Serre relations. 

The generators J ^ y can be written in several equivalent bases. The three most 
useful ones are the adjoint, bi-fundamental and auxiliary site bases. Let us present and 
compare them in what follows. 

3.2 Adjoint Basis 

In the adjoint basis the Lie and Yangian generators are denoted by and 3^^, re- 
spectively, where >1 = 1, . . . , A^^ for g = 0^-^)- The indices A transform in the adjoint 
representation of 

[j^,j;^]=F^^:^^. (3.1) 

The Serre relations are given by 

[3^^ [j^ y^]] + [y\ [j^, y^]] + [3^^, [j\ 3^^]] = J^, JO, (3.2) 

where {. . .} is the totally symmetric product (with unit weight) and 

A ABC _ 1 t^AGt^BH T^CJ (o o\ 

Here indices of the structure constant are raised and lowered by the (inverse) Cartan 
metric Cab, C^^. For the case of = 0^(2) the above cubic Serre relations are trivially 
satisfied and have to be supplemented by the quartic Serre relations 

[[3^^ 3^^], [J^, 3^^]] + [[3^^, 3^^], [J^ 3^^]] = {AiraF^^ + A%^,^Fr) {J^ J', y% 

(3.4) 

In all other cases the quartic relations follow from the cubic relations. 

We can use the latter property to unify the computation for all q\.{N) with N = 2 
and N > 2: In our case the proof of the cubic relation will make no reference to N 
whatsoever, explicitly or implicitly. Neither would a similar direct proof of the quartic 
relation make reference to N. Therefore, if the quartic relation holds for any one value 
of N, it will also hold for N = 2. Now because the above property tells us that the 
quartic relation follows from the cubic one for any > 2, it is sufficient to prove the 
cubic relation (for all N or for any N > 2). 

3.3 Bi-Fundamental Basis 

In the bi-fundamental basis for q = qI{N) the generators are denoted by J'^, y^, where 
a, 6 = 1,...,A'^. The adjoint transformation rule reads here 

{j^.y''A = ^iyi-^iyi (3.5) 

The commutator of two Yangian generators can be written explicitly, it yields the genera- 
tor y at the next level of the Yangian algebra phis a cubic combination of Lie generators 

[yh = ^imi - + Juu^ - j:JdJb- (3-6) 
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Likewise the higher-level generators (along with polynomials in the lower generators) can 
be obtained from repeated commutators of y. In the cubic Serre relations the y' terms 
drop out 

[Jb, y}]] - [yt, \ji yt\ = Ut, jfj^J! - J^jpi] 

+ ^){J2JtJ! - JPU^)- (3.7) 

The quartic Serre relation takes a similar form which we will not need explicitly here. 

To convert between adjoint and bi-fundamental basis one employs some tensors (tyi)c 
and (tc)"^ such that 

Jl = {tAtJ^, = {t^)lJ!, (3.8) 

and similarly for y. The two tensors are inversely related by (tA)c('^'^)fe = 

3.4 Auxiliary Site Basis 

The space of N x N matrices equals the space of qI{N) generators which equals the space 
of linear operators acting on an A^- dimensional vector space. Therefore we can use the 
latter as a basis for g[(A^)-adjoint generators. More explicitly we define 

J. = {J.)^CabJ'' = {JccfaJ^. (3.9) 

Here (Jlx)^ and (J^xTa qK^) generators acting on an auxiliary fundamental site 
labeled x. Consequently JJx is a linear operator that acts on the site x as well as on the 
spin chain (via J'^ or J^^). The adjoint transformation rules read 

\\Jxi\Jy\ \Pxyi3^x\i 

[Jx.yy] = [VxyM (3.10) 

where Vxy is the permutation operator. The cubic Serre relations are given by 

[Jx. [yy. yz]] - [yx. [Jy. yz]] = [X, Jz[Vyz. Jz\Jz\ + \Pyz. Jz[Vxz. Jz\Jz\ ■ (3.11) 

This basis has the advantage that one does not have to deal with explicit indices. 

3.5 Short-Range Chains 

In order to understand the structure of the Yangian for short-range chains, let us consider 
a specific known example: the Inozemtsev chain [T5l|16]. There, the generator y^ has 
the following simple structure 

y'' = J2l^i-k{X)FicJi'jf. (3.12) 

k<£ 

Here Cfc(A) are some coefficients whose precise expression will not be relevant. The 
coupling constant A is the coupling constant of the Inozemtsev model which is chosen 
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such that at A = the model reduces to the Heisenberg chain. For the Yangian generator 
this imphes Cfc(O) = 1. If one considers small values of A one finds that Cfc(A) = 1 + 0(A'^). 
The generator written as a series thus reads P^[^ 

oo 

_ 1 -pA q-B rrC , 1 \m -pA rrB q-C 

^ - 2^ 2^BC^k ^ 2^ 2^l-k^ ^BC^k J I ■ 

k<l k<i m=l-k 

oo m 

= X] l^Bc^kJf + X] X] XI I'^i^^^BcJ'k^J'k+e- (3-13) 

k<i m=l k 1=1 

In general we can write the above expression as 

oo 

y = yu + Y,^'y^''^^ (3-14) 

k=l 

where y^'^^ is a (local) operator which acts on no more than k + 1 adjacent sites of the 
chain and 3^bi is the bi-local operator 

iyu)t 

(3^bi)x 

We expect that the Yangian generator for general models of Sec. |2]is also of the form 
(13.141) . The bi-local piece 3^bi in fl3.15p is universal while the details of the model will be 
encoded in the precise expressions for the local pieces y^''\ Our expectation agrees with 
the generic asymptotic expression for the Yangian generator in perturbative short-range 
chains 

y''''J2,^jFBcJk''JF- (3.16) 

This expression explains why there should be no perturbative corrections to the bi-local 
piece J^bi in flS-lSp : the Lie symmetry generators are exact, they do not receive cor- 
rections in A, and therefore their bi-local product is asymptotically exact. In a chain 
with perturbatively modified Lie generators, such as [24j, the bi-local pieces also receive 
perturbative corrections, albeit given through the corrections of the Lie generators. The 
local terms y^''^ can be understood as a proper regularization of the above expression 
when the two Lie symmetry generators J'^ are close, i.e. within the range of the Hamil- 
tonian. These certainly depend on the details of the Hamiltonian if they are to commute 
with it. 

4 Construction of the Yangian 

In this section we will construct the Yangian generators for the short-range spin chain 
model introduced in Sec. [2l In their construction we will require that the Yangian 
generators commute with the local charges (including the Hamiltonian) and that they 



El pA rrB q-C 
k<l 2^BcJk Jl , 
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oi ... a„ 




1 



Figure 2: Graphical representation of a local adjoint interaction with one auxiliary site and n 
spin chain sites 

satisfy the Serre relations. The existence of this conserved Yangian shows that the 
associated Hamiltonian describes an integrable system and suffices to show perturbative 
integrability. In the construction of the local charges and the Yangian, we will assume a 
general form with many undetermined coefficients which are fixed by imposing the above 
constraints. 

4.1 Structure of the Perturbative Yangian 

The symbol y^'^^ contains the order A*^ perturbations which are local adjoint interactions 
of range k + 1. Note that we explicitly include contributions which act trivially on the 
auxiliary space (invariant interactions) as well as terms proportional to the Lie generator, 
J'. In principle we are free to drop these contributions: local invariant interactions 
which commute with the Hamiltonian define the local charges. Perturbing the Yangian 
generator by a local charge or by does not alter its adjoint transformation or the Serre 
relations. We nevertheless would like to keep these degrees of freedom explicit, because 
it simplifies the counting. 

This means that the Yangian generator can be written as some bare generator y plus 
some interactions proportional to the Lie generator or the conserved charges 

oo 

y = y + ly,j J + ly,c ^ + 5^ ly,s Qs- (4. 1) 

In order to match the perturbative interaction range we have to demand that 

ly,j = 0{X'), 73;,^ = 0(A°), 7y,, = 0(A^-i). (4.2) 

4.2 Notation 

As explained in Sec. 13.41 the Yangian is an adjoint operator which can be thought of 
as acting on the spin chain and an auxiliary site which is not on the spin chain. The 
notation introduced in Sec. 12.21 does not suffice to write these operators, so here we 
extend the notation to account for adjoint operators and bi-local operators. 

We denote an adjoint local operator which acts locally on n adjacent sites of the 
chain by the symbol [ CLq I I Cli , Cl2) • • • ) CLn 

]. A graphical representation of this operator is 
given in Fig. [21 Here the labels Oq, Oi, . . . , a„ are distinct integers from 1 to n and one 
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of them is the symbol 'x' representing the auxiliary site. The label ao denotes the site 
which is mapped to the auxiliary site and the other n labels describe which sites are 
mapped to the sites on the chain. As before there is an implicit sum over all n adjacent 
sites of the chain. For example an operator . . .] acts identically on the auxiliary site 
and therefore is equivalent to the invariant operator [...]. The operator would 
interchange a site on the spin chain with the auxiliary site. In fact this is just how the 
Qi{N) symmetry generator acts in (12. 4p and we have found a way to represent it in our 
notation 

J. = [l\\x]. (4.3) 
Simplifications as in fl2.7p also apply to adjoint operators 



ai, 



, a„,ra + ll 



ai, 02, 



[ao + l|| l,ai + l,a2 + l,...,a„ + l], (4.4) 



where in the last term the auxiliary label 'x' is clearly not shifted by 1. It can be shown 
that there are (n + 1)! — n! + 1 local adjoint operators of range n. This number coincides 
with the number of local invariant operators of range n + 1. 

To represent bi-local adjoint operators we use the symbol [ao| | ai, . . . , a^l a^+i, . . . , On]. 
A bi-local operator acts on k and n — k adjacent sites which can be at arbitrary distance 
on the spin chain. The first block of sites will always be left of the second block. Again 
ao denotes the site which is mapped to the auxiliary site. The labels ai, . . . , denote 
to which sites the first block of sites is mapped and Ofc+i, . . . , represents the second 
block. Again there are implicit sums over the positions of k and n — k adjacent sites such 
that the first block is always towards the left of the second block. In Fig. [3] we show a 
generic bi-local operator. In fact one may think of the bi-local adjoint operator as the 
following sum of local operators 



[aoll fli. 



1=0 
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''ki 



k + l,...,k + i,a[ 



(4.5) 



The primed labels equal a'j = aj + £ when aj > k and a'j = aj otherwise. The trouble 
with this expression is that it contains infinitely many terms and cannot be represented 
easily in a computer algebra. Instead one should implement directly the original bi- 
local interaction symbols. In this notation the bi-local piece of the Yangian fl3.15p is 
represented by 

(3^bi). = [l||2|a;]-[2||x|l]. 
The simplifications for bi-local operators are slightly more elaborate 



(4.6) 



oi, . . . , Ofcl afc+i, . . . , a„, n + 1] = [aol I fli, . . 
[ao + 1| 1 1, ai + 1, . . . , Ofc 1| cifc+i + 1, . . . , a„ + 1] = [aol I ai, . . 

[a'ol I a[, . . . , a'j^, k + 1\ a'j^^^^, . . . , 4] = [^o| I «i, • • 

- [floll Oi, 

[ao||a'i,...,a'fc|fc+l,a';,+i,...,a^] = [ao||ai,.. 

- [aoll ai. 



a„ 



a. 



nj ) 

an] 

■ ; Oiji] , 

an] (4.7) 
. , aJ. 



The latter two rules take into account that there is always at least one site in between 
the two blocks. The primed labels equal a^ = a^ + 1 if a^ > A; + 1 and a^ — CL£ otherwise. 
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ai ... Ofe cik+i ... a„ 




1 ... /c + n 



Figure 3: Graphical representation of a bi-local adjoint operator. 
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ansatz for Yangian 3^ 
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constraint from [H, 3^] = 
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15 
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undetermined coefficients 




2 
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26 


jy^s (mixing of Yangian) 
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5 


intrinsic parameters 
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5 


21 



Table 3: Number of degrees of freedom when starting with conserved local charges. 

4.3 Starting from Conserved Local Charges 

In the first approach we start with the results of Sec. [21 i.e. with a Hamiltonian Ti and a 
conserved charge Q. We then assume the most general expression for the Yangian fl3.14l) 
with (fc + 2)! — (fc + 1)! + 1 undetermined coefficients for the local adjoint operators at 
0{X^). Finally, we require that the Yangian generator commutes with the Hamiltonian 
and the conserved charge, as well as that it satisfies the Serre relations. This calculation 
is done explicitly to first order in A in App. [Bl at higher orders it proceeds in a very 
similar but not very illuminating fashion. In Tab-Owe list the number of free parameters 
and the constraints placed on them at each perturbative order. The resulting expression 
for the bare Yangian generator 3^ is presented at the end of the paper in Tab. [71 We 
discuss it further below. In App. [B] we construct this Yangian to C(A) as a concrete 
example. 

The first important comment is that the system of equations can be solved, thus 
we have found a perturbatively integrable system. The second outcome of this compu- 
tation is that the requirement imposes no further constraints on the coefficients of the 
Hamiltonian or of the local charge, but only on the coefficients of the Yangian generator. 
Therefore the Hamiltonian found in Sec. [2] is indeed integrable for the most general set of 
parameters. This confirms the observation [H[n] that for the present type of short-range 
spin chain models it is sufficient to impose the existence of merely one conserved local 
charge, Q, to obtain an integrable Hamiltonian. The surprising third comment is that all 
the constraints originate from the conservation of the Yangian generator; commutation 
with the conserved charge and the Serre relations follow automatically without further 
constraints on the coefficients. Finally, the Yangian is completely fixed up to the obvious 
deformations in fl4.ll) . 



14 





± 


A" 




A^ 


A^ 


ansatz for Hamiltonian Ti 




2 


5 


19 


97 
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73;, s (mixing of Yangian) 
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intrinsic parameters 
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21 



Table 4: Number of degrees of freedom when starting with conserved Yangian. Note the same 
number of intrinsic parameters as in Tab. [3l 

4.4 Starting from Conserved Yangian 

We have seen that merely one conserved local charge already appears to be sufficient to 
restrict to an integrable model. A useful question is whether different subsets of the full 
set of integrability conditions already guarantee integrability. In other words, we will 
now impose the various integrability conditions in a different order. 

Here we start with a conserved Yangian for an otherwise unconstrained Hamiltonian, 
i.e. we do not impose the existence of a conserved local charge Q. We use the same 
expansion for the Hamiltonian and for the Yangian as before. By imposing [7i, 3^] = 
we fix some coefficients of the Hamiltonian and the Yangian. We keep track of the 
number of undetermined coefficients at each perturbative order in Tab. HI From the 
total number coefficients in the Hamiltonian and Yangian we have to subtract the number 
of constraints and the number of degrees of freedom of choosing Ti and 3^ as a linear 
combination of all local suitable local charges. The remaining number of undetermined 
coefficients for the integrable system coincides nicely with the number obtained in Tab. [T] 
by taking all constraints into account. Thus, already demanding a conserved Yangian 
effectively constrains to an integrable system. Note that the Yangian and Hamiltonian 
here are identical to those of the previous section, this is not just a coincidence of 
numbers. 

We can now go on to fix the conserved charge by making it commute with the Yan- 
gian. This step does not introduce further constraints on the Yangian or the Hamiltonian 
obtained previously, it merely adjusts several coefficients of Q. This charge then auto- 
matically commutes with the Hamiltonian. 

4.5 Starting from Serre Relations 

Finally, we would like to see whether the Serre relations by themselves define the same 
integrable model. In fact, the Serre relations guarantee integrability, so we have to show 
that we can construct a suitable Hamiltonian which commutes with the Yangian. This 
procedure is quite interesting as it makes no reference to local charges at all. Conse- 
quently, the coefficients 7^,3 governing the linear combinations of local charges will enter 
only at a later stage. 

The Yangian is expanded as in (13.141) . Only 3^bi is specified, the perturbations given 
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21 



Table 5: Number of degrees of freedom when starting with Serre relations. Note the intrinsic 
parameters are identical to those in Tab. I3|4[ 



by adjoint operators are left unconstrained. We then require this form to satisfy the 
Serre relations in (13.111) to each order. Since receives no corrections the r.h.s. is zero 
at subleading order so we can use simplified relations which only have commutators of y 
and V, see (IA.2p . In Tab. [5] we keep track of the degrees of freedom. At each perturbative 
order we take the number of coefficients in the Yangian ansatz and subtract those fixed 
by the Serre relations. The undetermined degrees of freedom agree with those from the 
two other approaches above. 

The existence of a Yangian proves integrability if we can find a local operator that 
commutes with it; we need to construct a suitable Hamiltonian. We make the same 
ansatz for the Hamiltonian as before. In counting the degrees of freedom we must 
again subtract the constraints and the degrees of freedom used in choosing 7i as a 
linear combination of all suitable local charges. The constraint of [H, 3^] = completely 
specifies the Hamiltonian but gives no additional constraints on the Yangian, so we 
have the same number of undetermined coefficients as above. We can also construct a 
local charge Q in the same way and obtain similar results, i.e. no further constraints on 
the Yangian generator. This charge automatically commutes with the Hamiltonian as 
expected. 



4.6 Summary 

In the first method, we first constructed the Hamiltonian, 7i, and the conserved charge, 
Q, by requiring that they commute: [H, Q] = 0. We then constructed the Yangian gener- 
ator 3^ by demanding that it is conserved: [3^, H] = 0. We finally checked that 3^ satisfies 
the Serre relations and that it commutes with Q. In the second approach, we started 
with generic 3^, H and required [3^, H] = 0. We then showed that this Yangian generator 
satisfies the Serre relations. In the final approach we constructed a Yangian generator 
by requiring that it satisfies the Serre relations. We then constructed a Hamiltonian by 
requiring [3^, 7i] = 0. 

All three methods turn out to give the same number of undetermined parameters. 
Moreover the obtained Hamiltonian (Tab. [6]), conserved charge (Tab. [8]) and Yangian 
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(Tab. Ui> fully agree in the three approaches. This means that, at least for this system, 
the three approaches are equivalent. Curiously, imposing a single integrability property 
such as 

• the existence of one conserved charge, 

• the existence of a conserved Yangian-like operator, 

• a Yangian obeying the Serre relations, 

seems to imply full integrability, at least within the current class of models. In particular, 
we observe in the three cases that 

• many conserved charges [HE] and a conserved Yangian exist, 

• the Yangian obeys the Serre relations, 

• a local Hamiltonian which commutes with the Yangian exists, 

even without imposing further constraints on the existing operators. This is clearly not 
a rigorous proof, but it supports the claim in [HE] that for the present class of short- 
range spin chains, one additional conserved local charge is sufficient to prove that a 
Hamiltonian is integrable. 

5 Conclusions and Outlook 

In this paper we have proved the (third-order) perturbative integrability of the 0l(A^) 
long-range spin chain Hamiltonian proposed in [6j by means of Yangian symmetry. In- 
tegrability has long been assumed in the use of the asymptotic Bethe ansatz for such 
systems which predominantly come to use in the planar AdS/CFT correspondence [HI?]. 
This also supports the observation of [2S] that merely one conserved local charge is often 
sufficient to constrain a spin chain system to be integrable. 

We believe that a similar Yangian symmetry can be used for algebras other than 
0[(iV) and for different representations. Work has already been done to find the Yangian 
of a su(2|l) long-range spin chain [21] . Studying larger algebras involving fermions, gauge 
fields and derivatives (in the AdS/CFT context) leads to complications on the spin chain 
side such as fiuctuations of the length [27j for which new machinery will be required. 
In general it would be interesting to find out to what extent long-range integrable spin 
chains can be extended. 
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Figure 4: Graphical representation of the commutator of two local invariant interactions; here: 
[[2,1], [3, 1,2]]. 
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Figure 5: Pairwise cancellation of contributions to [J'x,Ly]. The cancellation extends to all 
terms when J'x interacts with a spin chain site a which maps to/from another spin chain site 
(i.e. those which are not labeled by 6 or y). 

A Simplifications of Commutators 

In the construction of the Yangian we come across many different commutators. In or- 
der to make the problem tractable we utilize several cancellations which simplify these 
commutators. To see these cancellations it is useful to think of the operators diagram- 
matically as in Fig. [T] and Fig. O We will describe the first simplification in detail, the 
others can be done in a similar fashion. 

When taking commutators we must consider the placements of the operator on the 
spin chain. As mentioned above there are implicit sums for these operators. When we 
take a commutator we need only consider the terms where the operators overlap as in 



Commutators of Local Adjoint Operators. The first commutator we will consider 
is [J'x, Ly] where Ly is an arbitrary adjoint local operator acting on the spin chain and 
on the auxiliary site y. The symmetry generator for qI{N) is defined in (14.31) . We should 
think of both J'x and Ly diagrammatically as in Fig. [2l albeit J'x acts on only one site 
on the spin chain. When we commute Jx with Ly there are two possible contributions: 
either J^ interacts with a site on Ly which does not map to/from the auxiliary site as 
in Fig. [5] or with a site which does map to/from the auxiliary site as in Fig. [6l We see 
there is no contribution from the former but there is a contribution from the latter. 

We can also consider [Pxy,Lx], as in Fig. [71 and we see it is equivalent to [Jx,Ly]. 
Hence we have the relation 



Note that if we substitute Ly with Jy or yy in ( lA.ll) we get the adjoint transformation 
rules which justifies that we call Ly an adjoint operator. 



Fig.H 
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X y 




b h 
Figure 6: The remaining contribution to [JxiLy]. 




b b 
Figure 7: The commutator [Vxy,Lx]. We see that this is identical to Fig. [6l 

Commutators Involving Yangian Generators. The simplification above is a sim- 
ple one but the same technique can be used on more complicated commutators involving 
bi-local Yangian generators in fl3.14p . In the construction of the Yangian generator we 
need to perform commutators of the sort [(3^bi)x, L] where L is now a local invariant 
operator. When only one leg of (3^bi)x interacts with L we get a cancellation similar to 
the one in Fig. [51 This is also clear from fl3.15p since the commutator is equivalent to 
[J'x,L,] which is zero. Thus there is only a contribution when both legs of interact 
with L. This simplifies our calculations since the outcome of [(3^bi)a;5 L] will always be a 
local adjoint operator (which is simpler than a bi-local adjoint operator). 



Serre Relations. To implement the cubic Serre relation (13.111) it is helpful to first 
simplify it for our purposes. We use the adjoint transformation (I3.10p to convert some 
generators to permutations of auxiliary sites 

-[v.., [yy,y.]] - [Vxy, [y.,y.]] - [Vy., [y.,yy]] 

[j7x5 JT^ [T^yX, JT^;] jT^j ~l~ ["^yz ) jT^ [T^^Z ? J^z] jT^ j . (A. 2) 

The permutations are somewhat simpler because they do not involve sums over sites of 
the spin chain. 

We can now expand (1A.2P using the expansion of the Yangian generators in (13.140 . At 
leading order we see that flA.2p involves commutators of two bi-local generators. However 



it is clear that the equation holds at leading order because this is just the undeformed 
Yangian generator of [1]. At subleading orders the r.h.s. of ( ]A.2p is zero since for 
g[(A^) receives no perturbative corrections. We see that the subleading corrections of 
the l.h.s. of ( ]A.2p will include commutators of the form [(3^bi)a;, 3^^'^''] and [yx\yy''~^^]. 
Terms of the first type could give a bi-local contribution but fortunately all terms when 
only one leg of the bi-local interacts with the adjoint operator sum to zero in flA.2p . 
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In fact, this result holds for the commutator of 3^bi with any adjoint operator, not just 
those corresponding to deformations of the Yangian. Thus there are contributions from 
[{yhi)x,yy'^] when both legs of the bi-local interact with the adjoint operator and from 
[yi^\yi'''^^] which both give local contributions. 

B Sample Calculation at Next-to-Leading Order 

Here we shall perform explicitly to first order in A the calculation which is described in 
Sec. M and Sec. 

Local Charges. We start with the most general Hamiltonian, Ti, and conserved 
charge, Q. At leading order 7i has length 2 and Q has length 3 

H(°) = ao,i[l]+ao,2[2,l], 

= 6o,i[l] + \2[2, 1] + &o,3[2, 3, 1] + 6o,4[3, 1, 2] + bo,,[3, 2, 1], (B.l) 

where aj^k is the coefficient of the k^^ term of Ti.^^^ and likewise bj^k for Q^''\ We now 
impose the constraint ['H^'^\ Q^'^^] = 0. This commutator is 

[n^'\ = ao,2(&o,3 + &o,4) ([2, 4, 1, 3] - [3, 1, 4, 2]) 

+ ao,26o,5 ([2, 4, 3, 1] + [4, 2, 1, 3] - [3, 2, 4, 1] - [4, 1, 3, 2]) (B.2) 

The constraints on the coefficients read 

bo,3 = -&o,4, &o,5 = 0, (B.3) 

which leaves us with 

= ao,i[l] + ao,2[2,l], 

= + &o,2[2, 1] + &o,3[2, 3, 1] - 6o,3[3, 1, 2]. (B.4) 

Now let us count the degrees of freedom. There are 2 unique local operators of length 2 
and 5 of length 3. There are 2 constraints from imposing [7^(0), Q(o)] = 0. According to 
( 12.101) . we can rewrite Ti.^^^ and Q^^^ using linear combinations of the bare charges. This 
allows us to fully specify Qf^ and Qf^: 

7^(o) = (ao,i + ao,2)/:-ao,2Qf, 

Q^°^ = (Vi + bo,2)C - 6o,2Qf + 2z6o,3Qf (B.5) 

with 

^ = [1], 

= [i]-[2,i], 

= -f([2,3,l]-[3,l,2]). (B.6) 
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To 0{\) we increase the length of each operator by one. H^^^ has 5 possible terms and 
Q^^^ has 19 possible terms. We now require H and Q commute to 0{X) which imposes 
the constraint 

[7^(0)^ g(i)] + [7^«,Q(0)] =0 (B.7) 

This imposes 16 constraints and gives 

H« = ai,i[l] + ai,2[2, 1] + ai,3[2, 3, 1] - ai,3[3, 1, 2] + ai,5[3, 2, 1], 

= bi,i[l] + 6i,2[2, 1] + 6i,3[2, 3, 1] - 6i,3[3, 1, 2] + (6i,io + ai,5&o,2/ao,2)[3, 2, 1] 

- 6i,io[2, 3, 4, 1] + 6i,io[2, 4, 1, 3] + 6i,io[3, 1, 4, 2] - 6i,io[4, 1, 2, 3] 
+ (ai,5V3/ao,2) [2, 4, 3, 1] + (01,560,3/^0,2) [3, 2, 4, 1] 

- (ai,56o,3/«o,2)[4,l,3,2] - (01,560,3/00,2) [4, 2, 1, 3]. (B.8) 

As before we can rewrite these charges in terms of linear combinations (l2.1Up of the bare 
charges 

7^(1) = ^ + ai,2 + ai,5)/: + (-ai,2 - 4ai,5)Qf - ao,2Q5'^ + 2zai,3Q^°\ 
Q^^^ = (&1,1 + &1,2 + 61,10 + ^0,201,5/00,2)'^ 

+ (-61,2 - 261,10 - 460,201,5/00,2) - 60,2 Q?^ 

+ (2^61,3 + 12^01,560,3/00,2) Qf + 226o,3Q^'^ - 36i,ioQf , (B.9) 

where 

= (ai,5/ao,2)(-3[l] +4[2,1] - [3,2,1]), (B.IO) 

Qs'^ = 1(^1,5/00,2) (6[2, 3, 1] - 6[3, 1, 2] + [4, 1, 3, 2] + [4, 2, 1, 3] - [2, 4, 3, 1] - [3, 2, 4, 1]) , 

2f = I (-[1] + 2[2, 1] - [3, 2, 1] + [2, 3, 4, 1] + [4, 1, 2, 3] - [2, 4, 1, 3] -[3,1, 4, 2]) . 

We see that there is one degree of freedom left in describing the bare charges: 01,5/00,2- 
This counting corresponds well with Tab. [H The expressions flB.6IIB.10p agree with the 
expansion of Tab. [UHHl and we find an expression for the function ao{X) 

ao(A) = A(oi,5/oo,2) + 0(A2). (B.ll) 

Note that the results in Tab. and Tab. [H] use the coefficients described in Sec. 12. 4j 
these coefficients distinguish the physical degrees of freedom from the gauge degrees of 
freedom. 

For completeness we will give all the coefficients specified in fl2.10p to order A. Com- 
paring (IB.SP and ( IB. 91) to f l2.10p we get the following 

72, £(A) = Oo,l + Oo,2 + A(Oi,i + Oi,2 + 01,5) + C(A^), 

72,2(A) = -00,2 + A(-oi,2 - 401,5) + C(A2), 
72,3(A) = A(2mi,3) + 0(A2), 

73, £(A) = 60,1 + 60,2 + A(6i,i + 61,2 + 61,10 + ^0,201,5/00,2) + C'(A^), 
73,2(A) = -60,2 + A(-6i,2 - 261,10 - 460,201,5/00,2) + C(A^), 
73,3(A) = 2i6o,3 + A(2i6i,3 + 12201,560,2/00,2) + C(A^), 

73,4(A) = A(-36i,io) + 0(A2). (B.12) 
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Yangian Generator. Now that we have Ti, and Q to first order we can construct the 
Yangian, 3^. To leading order this is given by the bi-local operator plus invariant terms 
of length 2 which gives the Yangian ansatz 2 degrees of freedom. 

= [111 2| x] - [2\\x\ 1] + co,i[a;|| 1] + co,2[l||x]. (B.13) 

We then consider the constraint 

[{y^'^).,n^'^=o, (B.14) 

which is satisfied automatically. According to (14. ip we can rewrite 3^^"'' as a linear 
combination of the bare Yangian, the bare charges and the symmetry generator 

(3^^°^). = (3^^°^). + co,iC + co,2J., (B.15) 

where 

(3^(°)). = [l||2|a:]-[2||x|l]. (B.16) 

Thus 3^*^"'' is specified exactly and there are no new constraints on H.^^^ 
To 0{X) we have the following most general ansatz 

{y^'^)x = ci,i[x\\l] + ci,2[l|| x] + ci,3[a;|| 2, 1] + ci,4[l|| 2,x] + ci,5[2|| a;, 1]. (B.17) 

Requiring Ti. and 3^ to commute to 0{X) is equivalent to the condition 

[{y^%,n^'^] + [(3^('))x.,H(°)] = (ao,2Ci,4-2ai,5)([l||3,x,2] - [2||3,l,x]) 

+ K2Ci,5 + 2ai,5)([3|| 2,x, 1] - [2|| x,3, 1]) 
= 0, (B.18) 

which yields the constraints 

ci,4 = 2ai,5/ao,2, Ci,5 = -201,5/00,2- (B.19) 
Using (14.11) we now rewrite 

37(1) 

as a linear combination of the bare Yangian, bare charges 

and symmetry generator 

(3^^^^).. = (3^^^^).. + (Cl,l + ci,3)/: - Ci,3Qf + Ci,2 J., (B.20) 

with the next-to-leading order bare Yangian 

(3^«), = 2(ai,5/ao,2)([l||2,a;] - [2||a;,l]). (B.21) 

As with the leading-order case, 3^'-^^ is specified exactly and there are no new constraints 
on H^^^ This counting is continued to 0{X^) in Tab. [31 

For completeness we once again give the coefficients specified in (14. ip . Comparing 
(IRTsD and (lOoD to (gl]) we get 

7y,j(A) = co,2 + A(ci,2) + 0(A'), 
ly,cW = co,i + A(ci,i + ci,3) + 0{X^), 

7y,2(A) = A(-ci,3) + 0(A2). (B.22) 
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Serre Relations. Finally we have to show that this Yangian satisfies the cubic Serre 
relation (lA.2p . The details of this are not very illuminating so they are omitted. This 
Yangian we have constructed also commutes with Q automatically. So given the system 
formed by requiring one conserved local charge [H, Q] = 0, we can construct a Yangian 
which imposes no further constraints on the system but guarantees integrability. 

We were only explicit here up to order A but the higher orders follow in the same 
way. The other methods are also done in a similar fashion, albeit with the constraints 
imposed in a different order. 
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Q2(A) = ([1]-[2,1]) 



+ ao(A)(-3 [1]+ 4 [2,1] -[3, 2,1]) 

+ ao(A)'(20 [1] - 29 [2, 1] + 10 [3, 2, 1] - [4, 2, 3, 1] 

+ [3, 1, 4, 2] + [2, 4, 1, 3] - [4, 1, 2, 3] - [2, 3, 4, 1]) 
+ f ai(A) (6 [2, 3, 1] - 6 [3, 1, 2] + [4, 1, 3, 2] + [4, 2, 1, 3] - [2, 4, 3, 1] - [3, 2, 4, 1]) 
+ (-4 [1] + 8 [2, 1] - 2 [3, 1, 2] - 2 [2, 3, 1] 

-2[2,l,4,3]-2[4,l,2,3]-2[2,3,4, 1] 

+ 2[3,l,4,2] + 2[2, 4,1, 3]-2[3, 4,1,2] 

+ [4, 1, 3, 2] + [2, 4, 3, 1] + [4, 2, 1, 3] + [3, 2, 4, 1]) 
+ ze+,(A)([2,4,l,3]-[3,l,4,2]) 

+ ze+2(A) ([4, 2, 1, 3] + [2, 4, 3, 1] - [4, 1, 3, 2] - [3, 2, 4, 1]) 

+ e>(A') 

Table 6: Normalized Hamiltonian printed up to second order. 



5>. = ([l||2|x]-[2||x|l]) 

+ 2ao(A)([l||2,,T]-[2||,T, 1]) 

+ 2ao{\f (2 [2| I x, 1] - 2 [1 1 1 2, x] + [1 1 1 3, 2, x] - [3| | x, 2, 1]) 
+ iai(A) (2 [l||2,3,a;]+2[3||a;,l,2] 

-[l||3,x,2]-[2||x,3,l]-[3||2,x, l]-[2||3,l,x]) 
+ /?2,3(A) (2 [l||2,x]-2[2||a:, 1] 

- [1 1 1 3, X, 2] + [2| I X, 3, 1] + [3| I 2, X, 1] - [2| I 3, 1, x]) 
+ 2ze+i(A)([3||2,l,a:]-[l||x,3,2]) 

+ 2ie+2(A) ([2| I 3, 1, x] + [3| I 2, X, 1] - [1| I 3, x, 2] - [2| | x, 3, 1]) 
Table 7: Normalized Yangian generator printed up to second order. 
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(A) = |([3,l,2]-[2,3,1]) 

+ |ao(A) (6 [2, 3, 1] - 6 [3, 1, 2] + [4, 1, 3, 2] + [4, 2, 1, 3] - [2, 4, 3, 1] - [3, 2, 4, 

+ |q;o(A)'(46 [3, 1,2] -46 [2,3,1] 

+ 12 [3, 2, 4, 1] + 12 [2, 4, 3, 1] - 12 [4, 2, 1, 3] - 12 [4, 1, 3, 2] 
+ 2 [2, 4, 1,5, 3] + 2 [2, 3, 5, 1,4] + 2 [3, 1,4, 5,2] + 2 [5, 1,2, 3, 4] 
-2 [4, 1,2, 5, 3] -2 [3, 1,5,2, 4] -2 [2, 5, 1,3, 4] -2 [2, 3, 4, 5,1] 
+ [5, 2, 3, 1,4] + [5, 1,3, 4, 2] + [5, 2, 1,4, 3] 
-[2, 5, 3, 4,1] -[4, 2, 3, 5,1] -[3, 2, 5, 4,1]) 

+ iai(A) (-20 [1] + 24 [2, 1] - 8 [4, 1, 2, 3] - 8 [2, 3, 4, 1] - 4 [4, 2, 3, 1] 
+ 6 [3, 1,4,2] + 6 [2, 4, 1,3] + 2 [4, 3, 1,2] + 2 [3,4,2,1] 
+ 2 [5, 1, 3, 2, 4] + 2 [2, 4, 3, 5, 1] - 2 [2, 5, 3, 1, 4] - 2 [4, 1, 3, 5, 2] 
+ [5, 2, 1, 3, 4] + [3, 2, 4, 5, 1] + [5, 1, 2, 4, 3] + [2, 3, 5, 4, 1] 

- [3, 1, 5, 4, 2] - [4, 2, 1, 5, 3] - [2, 5, 1, 4, 3] - [3, 2, 5, 1, 4]) 

+ i/92,3(A) (-4 [3, 1, 2] + 4 [2, 3, 1] + 4 [4, 1, 2, 3] - 4 [2, 3, 4, 1] - 2 [4, 2, 1, 3] 
-2 [4, 1,3,2] -2 [3,4,2,1] + 2 [4,3, 1,2] + 2 [3,2,4, 1] + 2 [2, 4, 3, 1] 
-4 [4, 1,2, 5, 3] -4 [3, 1,5,2, 4] -4 [2, 3, 4, 5,1] -4 [2, 5, 1,3, 4] 
+ 4 [2, 4, 1, 5, 3] + 4 [2, 3, 5, 1, 4] + 4 [3, 1, 4, 5, 2] + 4 [5, 1, 2, 3, 4] 

- [5, 2, 1, 3, 4] - [3, 1, 5, 4, 2] - [5, 1, 2, 4, 3] - [3, 2, 5, 1, 4] 
+ [2, 3, 5, 4, 1] + [4, 2, 1, 5, 3] + [3, 2, 4, 5, 1] + [2, 5, 1, 4, 3] 

+ 2 [3, 1,2, 5,4] + 2 [2, 1,5, 3,4] -2 [2, 1,4, 5, 3] -2 [2, 3, 1,5,4] 

+ 2 [4, 1,5,2, 3] + 2 [3, 5, 1,2,4] -2 [5, 1,3, 2,4] -2 [3, 4, 1,5,2] 

+ 2 [2, 4, 3, 5,1] -2 [2, 4, 5, 1,3]) 
+ |e+i(A) ([4, 2, 1, 3] + [3, 2, 4, 1] - [4, 1, 3, 2] - [2, 4, 3, 1] 

+ [3, 1, 4, 5, 2] + [2, 5, 1, 3, 4] - [4, 1, 2, 5, 3] - [2, 3, 5, 1, 4]) 
+ |<2(A) ([5, 2, 1, 3, 4] + [3, 1, 5, 4, 2] + [3, 2, 4, 5, 1] + [2, 5, 1, 4, 3] 

- [4, 2, 1, 5, 3] - [5, 1, 2, 4, 3] - [3, 2, 5, 1, 4] - [2, 3, 5, 4, 1]) 

+ 0(A^) 

Table 8: Normalized conserved charge printed up to second order. 
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